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Section 1 


Let us assume we must solve the vector equation 


x=Ax-*-y, (1) 

where y is the known; x is the unknown n-dimensional vector; A — 
n-dimensional matrix. TOie method of successive approximations in 
solving equation Cl) consists of the fact that we start with the 

arbitrary vector x° and then subsequently formulate the system of 

, iaAfn „ 1 2 k, . . . , where 

vectors x , x , . . . , x * * 

jf=Ax*~ 1 -*-y. (2) 

k 

If the sequence x converges to the limiting vector x, then 
x is the solution of equation (.1). However, the sequence con- 
vergence may be very slow, and the problem arises of extrapolating 

1 k 

the approximations already found of x , ...» x to the root x, so 
as to approximate it more closely, avoiding the repeated substitu- 
tion of (2). 

Let us examine the simplest case when A is the symmetric matrix. 
Let us assume , • • • , are the eigenvalues of A, and 

jX„ U| 9 u t , is the orthonormal system corresponding to 

it of eigenvectors. Subtracting the equation (1) obtained from 
(2) , we have 
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( 3 ) 




f—x-Atf-'-x). 


If 2 e <“<» » then 


tel 


tel 


tel 

and, in general, 

N 

jc i ‘— x= V| V e ,“#* 




(*> 


Repeating the estimates used when determining the eigenvalues 
numerically using the Mises method, we find: if |>,|=|V = ...:= 

W=x,IVil<i at 3 > 0, then 0 at and 

*?—x= X 1 f 2 *, c ( u, c,+j (^p) ~ X* 21 €* c. «„ K I — 1 

\ tel tel / tel 


Z§0 


(If one of the numbers c, , c«,..., does not equal 0). For 
large k, the expression x - x may be approximately represented 
by the sum of the eigenvectors corresponding to the eigenvalues 
which are largest in absolute value. In the case of a symmetric 
A and, consequently, real eigenvalues, one or two (opposite in 
sign) eigenvalues (of any multiplicity) may be the largest in ab- 
solute magnitude. The criterion for the convergence of the method 
described of sequential approximations is, as is known, the fact 
that the eigenvalue of the matrix A which is largest in terms of 
absolute value is less than unity. 


For purposes of simplicity, let us examine the case when there 
is a singular (of arbitrary multiplicity) eigenvalue of X which is 
the largest. Then 


^ — x 


X*tt, 


(5) 


where u is a certain eigenvector corresponding to the eigenvalue 
of X. Similarly, 


**♦' — u, 

x t + a —x~y. i +*u. 


( 6 ) 
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Thus 


**— X)tf, 

The vectors x*— x** 1 and jh- for large k are almost pro- 
portional, and the ratio of their components approximately equals 


1 — x* 


l+x. 


Thus, If the sequential approximations x *, J* 1 , are known, 
then we may approximately find X , and since we have the following 
from C5) and C7) 


If 

by determining X , we may approximately find x - x, and this also 
means x. 

If x Is very close to unity, then It is advantageous to use 
similar equations instead of C.6] 

jf+v —x — X' +i '' u, ^ 


from which we have /51 

— ).*»■) u, 

and, consequently, the ratio of the components x'—x** 1 ’ and 
x***'' determines 1 -+ 0 / , and this means X also. 

In addition. 


x* —x ~ — x>+ 0- 
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Let us now assume operator A has two of the largest eigen- 
values for X and *»=—>•. . Then 

J-x^x'iu+i-ifvj (8) 

Cu and u^ are the eigenvectors referred to the eigenvalues X and 
*^1, and In addition 

* l+ * — x — *' * +5 [u •+- (— 1 )‘ «»,], 

-.-<-! )>], (9) 


Therefore, the vectors **— jf +i and have approximately 

proportional components, and their ratios approximately equal 


from which we may determine X and in view of (8) and (9) find x 
approximately 


* r=Ti (** — *?*') ~ rzrr« (** - **+*)• 


Just as In the previous case. Instead of x \ J+* , we 

may use *• x*"- 


Section 2 


By way of an example, let us give the numerical solution of 
plane boundary value problem of the Laplace equation by the grid 
method. Let us first make some preliminary comments. 

Let us designate the following operators by 7* and D ■ D* 1 
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V* [*(*-♦- A, y) -*-*(*-- A, y)-+- 
-»-*(*, y + A) -»-*(*, y — A) — 4* (jr, y)J, 

Dz = D* *=qp 7**-* A, y)-f*(* — A, y) 

-*-z(x, y-t-h)-t z(x t y — A)], 


The operator V n is the finite-difference approximation /52 

of the operator V. Let us examine a plane square grid with one 

side of the square having the length h and the finite plane 

region Q, consisting of the squares of the grid. We shall give 

the function z for the grid only at the grid corners . Then 

V* 1 and D* 1 are changed into finite-dimensional operators (matrices). 

If the values of z at the corner of the grid (ih, kh), where i,k 

are whole numbers, can be designated by z. v , then the operator 

D n replaces z* » by the value 4- (:>* -*■*<.*-«). . Just 
1 » K h 

like the operator v, the operator V has only negative eigen- 
values “ •>—!/>••• . Th e operator D* 1 has the 

eigenvalue 

* , a* • . 

ft ® 1 — T '• ’ (10) 


while all ^*<1. 


If the number of Internal grid corners at Q equal N and all 
the functions u, subjected to the operation D h , have the values 
0 at the boundary corners of the grid <$, then V* 1 and D* 1 have N 

difference eigenvalues — — V. • . •• — i*i*. :*** with the 

general eigen function uj ; «*,•••, The eigen functions u^, u 2 > 

. which corresponds to the first eigenvalues — * : \ — V» •• approxi- 
mate the first eigen functions of the Laplace operator v with zero 
boundary conditions. 
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We should note that the eigen function u of the same operator 
with an eigenvalue which Is Inverse In sign corresponds to each 
of the eigen functions u of the operator D* 1 , where 

Cll) 


In the same way, If Dt then 


and since 


T' u 


*M. * 


<— I. t 


!♦. 




( 12 ) 


differ from 


!.»• U »-i. i * 




, i t °j_;. k» 

by the factor (— !)*♦*«(— 1)*~ : , we have the following from (12) 


7 [“<+!.*' 




t-1. 


t 


i.e., * 7116 functlon u is the eigen function of the 

eigenvalue — y. Therefore, the eigenvalue — y of this same opera- 
tor with the inverse sign corresponds to each eigenvalue y of the 
operator D h . The smallest eigenvalue ftr* corresponds to t’ie 
largest eigenvalue o.‘ , and corresponds to the 

eigenvalue . We have 


1 > ,*i* >*/>•••> = > ~ " ' ' > l| 


So that from all |.**j <i , we have the well-known fact that the /53 
Iteration process converges by means of the operator D ■ D h . 


We should note that the eigen functions u^, Uj,*.., corre- 
sponding to the first eigenvalues y t \ , approximate the first 
eigen functions of the operator V. The eigen functions u 2 , 

..., determined by means of formula (11), which corresponds to the 
negative eigenvalues — /*/» — •• » which are largest in absolute 
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value, are the singularities of only the grid approximation of 
the Laplace operator, but not the operator Itself. We shall oall 
these eigen functions parasitic. Their existence must be con- 
sidered when studying the convergence of certain approximation 
processes. 


Section 3 

The approximate solution of the boundary value problem for 
the domain of the Laplace equation by the grid method consists of 
finding u at the grid corners which satisfies the following equa- 
tion at the Internal comers 


and certain boundary value conditions at the boundary comers, for 
example, finding the values of u at these comers In the case of 
the Dlrlchlet problems. For simplicity, let us discuss this prob- 
lem. Starting with the arbitrary function w at the internal comers 

of the grid, which uses the given values at the boundary comers, 

12 k 

let us formulate the sequence of functions w , w ,..., w ,..., 
which satisfy the same conditions at the boundary corners, and 
the internal decisive conditions 


v 

The difference w — w satisfies the zero boundary value conditions. 
If u and u are the eigen functions of the operator D, corresponding 
to the largest and the smallest eigenvalue of \i and — y , then 

a>* — w—cu-i-eu-*-v, (l^j) 

where v Is the sum of the eigen functions referred to other eigen- 
values (if w° is selected as a smooth function, then w° — w Is a 
smooth function and the coefficient c, the coefficient with the 
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parasitic function u which changes sign for each pair of adjacent 
corners, is small). Repeating the previous discussion, we find 

w* — w — u* [cm 1)* is I 

- a 4 * I«# + (- i/m)* 

w* 44 — » ~,u*' w [«.-v-(-l)‘cal, 
w‘— - /(l - X*) [cu + (- lfci], 

-,*‘(1 - v)M • <— 

The functions u/—u/' 4 \ »*— «.*■*■* are approximately proportional, 
and their ratio approximately equals 






It is best of all to set 


Is- 


2w./--‘V) 

2(-*. ’ 


(16) 


where the sums are given for all internal corners of the domain Q. 
Determining A, we may use formula (15), from which we may find the 
approximate value of the function w. 


Actually, the calculations are as follows: We find 

w* 44 ; by the successive use of the operation D. We de- 
termine the sums 2“&*» for a11 internal corners, and 

in addition, we find 1 \ A using formula (16), we determine A us- 
ing formula (15) and find w. 


Let us examine the example given in the n Handbook for Numerical 
Solutions of Partial Differential Equations,” D. Yu. Panov (GTI, 
19**3, pp* 50-58) for the solution using the method of successive 
approximations of the equation Dw ■ w. 


Tables 1, 2 and 3 give the results of the 7, 9 and 11th itera- 
tions of the operation D. The values at each comer are located 
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: to the Toft and abovo It . The boundary values an only given in 
Table 1. The values at the comers which are symmetric with re- 

speot to the diagonal equal They are given only on 

one side of the diagonal («£/ designates the value at the comer 
of the intersection of the ith low and the jth column after the 
kth Iteration. Tables 1, 2, 3 give the values of •L* "O' 

Tables 4 and 5 give the differences and 

at all comers of the grid. We have 

V (o£ t — ««) K 2 (0.12 0.20 0.22 ■+- 0.17 0.21 0.36 + 

-f 0 39 0.27 0.43 0.24) -+■ 0.1 033 0.48 0.41 -«-0.14— 6.68, 

<^ = 2 (a07n-ai2 -«-0.12 0.10 H 0.12-1-0.21 -»- 0.22 
0.16 035 0.14) -4- 0.06 +■ ai8 -+• 0.27 ■+■ 0.24 0.08 = 3.95. 

Using the formula (16) 

l+>.'-{g-1.7J5, 

from which we have 


1 — >.* = 2 — 1.735 = 0365, T ~j 1 = 0 ^ j ==3.77. 

Table 6 gives the values of t =w l,~ 3>77 (*?. / "-"E t) [see 

formula (15)]. Table 7 gives the results of applying the averag- 
ing operation D to Table 6. Table 8 contains errors in hundreths /56 
of unity for the values of Table 8. This accuracy is achieved in 
the book of D. Panov after 26 iterations. 

For comparison, in Table 9 we give a table of errors in the 
values of w* after 26 iterations D according to the "Handbook.” 

Interpolation used l*i iterations. If w» use the formula 
w >., =w l>— (»’.>—<*,)> ln ^5), then the errors are somewhat 

smaller. 
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Section 4 


The equation Cl 3) nay be replaced by the equivalent equation 

(13*) 

where 


/V *>*■ »"*•*»). 

The eigenvalue! i*/ of the operator D^are connected with 
the eigenvalues p** of the operator D^ 1 by the formula y,*— 

Since all Jp|<l , then at *>0 all Ip.K* . The eueoeeelve use 
of the operator D a h gives a converging process, and to solve the 
system (13) or (13* ), we may use the successive use of the opera- 
tor D 0 h . 


If at— ><0, 


then 




N "TT 


However, since there may be values which rre very close to —1 /57 

among the eigenvalues i^ 1 , the eigenvalues p a h may be less than 
— 1. If — Is the first eigenvalue of the Laplace operator 
for our domain, for a small h It Is close to the first eigenvalue 
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h * , . 

— for the sane domain of the operator 7 i — V"*-*'! . 

In view of (10), we have 

~ 1 - -£ >|{ J»jr* «f*i* ^ v *, ~ 1 

which Is very close to — 1. The smallest of the Is 

, , A*. 

-j-/, 

rr* 

For a very small h and a fixed **»—£, , the smallest 

eigenvalue becomes less than —1. But then equation (**3) 
cannot be solved by Iteration of the operator D **, since It dis- 
turbs the necessary convergence condition of the iteration proc- 
ess (but not all of the eigenvalues of D Q h are less than 1 In 
absolute value). In this case, the convergence is called the 
parasitic eigen function. Only the operators D a h with positive 
a can compete with the operator D* 5 . 

The operation D ■ D* 1 consists of using the processes of 
multiplication Individually. In terms of simplicity in perform- 
ing the operation with the massive use of operations, operations 
which are reduced only to multiplication and division by tens are 
of Interest (l.e., moving the decimal). For example, this le the 
case with the operation for isj. 

7 (**.*♦! + *»♦».* ■*■**->.* 

The operation Is reduced to replacing the value of z ^ 

at each grid comer by the mean arithmetic value, and at 5 comers 

— by the given value and the 4 adjacent ones. The division of 
the nuriber by 5 reduces to doubling the number and dividing by 10, 
i.e., moving the decimal point. The successive use of the opera- 
tion # /t Is actually reduced to the successive use of the opera- 
tion of multiplication, which may be completely carried out on 
adding machines (from simple ones to tabulation machines). 
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If we do not take a square grid, but a rectangular grid. In 
which the sides are rectangles which are parallel to the z and y 
axes, referring both to j/S and 3% t then for the approximate solu- 
tion of the Laplace equation, we must find the function u at the 
comers of the rectangular grid which satisfies the following equa- 
tion 

where 


Df «*, ,= la i2 Ui+u ,+2^,-»-3o l>M . l + 3 

Oil k Is the value of u at the comer, which lies on the Inter- 
section of the ith vertical and the kth horizontal. Here the 
sequential use of the operation is only reduced to summation. 

Similarly, in the three-dimensional problem we may introduce /58 
the approximate solution of the boundary value Laplace problem by 
the grid method. One of the methods indicated above, or a combina- 
tion of them, may reduce all of the successive operations to sum- 
mation and division by 10, l.e., practically only to summation, 
and only to addition. For example, let us give this method. 

Let us consider the following three-dimensional grid: The 

plane xOy is divided Into triangles with equal sides with the side 
h with three series of parallel lines. The apexes of these tri- 
angles are corners of a plane grid. This grid is moved parallel 
along the z axis and is reproduced on the planes 

x—z H — nh^ (n= ± 1, — 2, • • •)• 

We obtain a three-dimensional grid, and the function u is given at 
the grid comers. Let us assume A is a grid corner Capex of the 
triangle) in the plane z ■ z R . 

Let us introduce the operation F* 1 , which changes the value of 
u(A) at the corner A, divided by 10, and the sum of the value u at 
6 comers of the same plane, which are adjacent to A, and the double 
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sum of the value u at 2 oorners lying above and below the point A 
in the plane z • z n+1 and z ■ z n _^. 

As may be readily seen, the approximate solution of the Laplaoe 
equation reduces to find the solution of the equation u » F* 1 u at 
the grid oorners. The solution of this equation by the method of 
successive approximation is reduoed to the successive use of the 
operation F h , i.e., only to summation. 

Let us assume we must find the three-dimensional solution of 
the Laplace equation with axial asymmetry. The axis of symmetry 
Is with the z axis, and we assume that the previous grid is con- 
structed so that one of the triangle apexes in the z ■ 0 plane 
(and in the parallel plane z ■ z R ) passes through the symmetry 
axis. 


The function u must assume equal values at the corners of 
each plane z • z u , which are equally removed from the axis. Since 
our grid has 6 symmetry planes passing through the axis, the number 
of oorners corresponding to the different values of u is approxi- 
mately 6 times less than the number of corners of the three-dimen- 
sional domain, and it approximates a plane grid using the number 
of unknown values of the three-dimensional problem thus solved. 


Section 5 


Let us examine the Sturm-Liouville equation 

l # (*) V (0] - P (x) y (.*•) - \y (x) ( 17 ) 

under the boundary value conditions, for example y(a)*0, y(bj*0. 
We will have R(x)> 0 at the segment ta» b). We shall also assume 
that on l«. *»1 P(x) ,'Q. 


Along with (1), let us consider the equation 




( 18 ) 


Let us assume u*(x) ,u 2 (x) , . . . , u k (x),... Is a sequence of eigen /59 
functions of the equation (17), corresponding to the values 

if P(*)<o on ta, b], then are nega- 
tive. If 

00 

?(*)=£**«*(*). , % 

*T. (19) 


then the solution of (18) with the initial value of w(x,0)«9.(x) 
is 


w(x,i)*= V c„ u (x)e >tl 




( 20 ) 


Since all X k < 0, then w(x,t) approximately strives to zero at 
t-+ c\>, and the terms c t u*(x)e v strive to zero more rapidly, the 
larger is k. If, for example, Cj f 0, then for a very large 
t the first term w(x,i)fcc t u' (x)e y '‘. will be the main term in (19). 




1 IT- 


U 


F 




- *-VW^1 


V77. 


Pig. 1 


Por a fixed rather large t, the function w(x,t) may become 
approximately proportional to the first eigen function 
with the relative accuracy which Increases as t increases. If 
we simulate the equation (18), then we may find the approximate 
values of the eigen functions of the equation (17), and the eigen- 
values . 


The equation (18) may be modeled for a linear electric circuit 
of resistances and capacitances. 
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p«r- 


jgsngegF? ^r^sf - •*•• srs ' i : " ‘ 



if 

s-f 


Dividing the segments la,b] into n equal parts with the length 

getting fCa+ih)-^, we replace equation (17) by 
its finite difference analog 

jr [R< (y i+ > - s.) - *.-« (9, - v.-i » ] - p * • ( 21 ) 

Along with (21), let us consider a system of ordinary equa- 
tions which approximate the partial differention equation (18) 

p- [£,(»,* (<) - w, (<)) + R .t (u.(<) - W )) ^ p < u,<<)= Cw t ‘ (<). ( 22 ) 

All of the previous considerations hold if we designate the 
’’eigen functions” of equation (21) by u k (k ■ 1, 2,.. ,). These 
functions are given only at the points a » lh, and we use X k to /60 
designate the corresponding eigenvalues which approximate the 
eigen functions and values of the initial equation (1) at h + 0. 

Let us formulate a linear scheme of resistances and capaci- 
tances with the corners k^ (i ■ 0, 1, 2,..., n) (Pig. 1). Each 
corner A^ is connected with the corners A,^ and A 1+1 by the con- 
ductivities Rj.j and R^, and currents with the conductivity Pj 
and the capacitance C emanate from it. The voltage equals 0 at 
the corners k Q and A n . The voltage at the corners A^(i * 1, 2, 

..., n —1) satisfies equation (22). At the corners A^, Ag,..** 

Aj, ...» A n _^ setting the arbitrary system of voltage, in view of 
the previous statements, after a certain time interval t, we ob- 
tain a system of voltages w, (t) which is approximately proportional 

1 1 

to the system of values for the first eigen function u at the 
points a + ih. If we give the voltage in hundreths of volts and 
measure them at the time t for a voltage of about one volt, then 
the measured voltages at the comers A^ give the approximate 
values of the first eigen function at the points a + ih. Thus, 
the sign and overtones — of other eigen functions, are retained. 

If the measured values of the voltage — when we multiply them, 
let us say, by 100 — are given at the corners A^ and are again 
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tis 


measured in the time interval t Q , then we may obtain the more 
"pure" value of the first eigen function (with weaker traces of 
the higher eigen functions). 

Let us assume we have obtained the first eigen function u, 1 
which assumes the value u^ at the points A^(i ■ 1, 2, . .., n — 1). 
Let us now give the system of voltages w at the corners A^, so that 

i.e., the initial distribution represents the function w which is 
orthogonal to the first eigen function u*. For example, let us 
assume that at all of the corners except two, w ■ 0, and at the 
corners and A^ 


— u t l : u, 1 . 

Similar measurements would have to give the second eigen 
2 

function u . The given function w is orthogonal to u', and 
therefore the coefficient c^ ■ 0 in the representation of (3) 
when <J> is replaced by w. Now, with an increase in t, the second 
term will decrease the most slowly 

e. u* (x) 1 

2 

with the following eigen function u . 

The coefficient c^ does not precisely equal zero for u', al- 
though it is small, but since the first term decreases more slowly 
than the others, the measured system of values at the time t does 
not give an expression which is approximately proportional to the 
second eigen function. Generally speaking, it may be regarded as 
a linear combination of the first two eigen functions. However, 
the first eigen system is known. Therefore, we may determine the /6 1 
second from this linear combination. 
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The transition to the third eigen function is similar. 



Let us assume />(*)^° in the Interval (a, b). If this 
condition is violated, then it is sufficient to replace the co- 
efficient P(x) in (17) by the coefficient 

so that it is positive in this Interval. This substitution leaves 
the eigen functions of equation (17) unchanged, and decreases all 
of its eigenvalues by the nuriber C. 



We may use a similar method to model the eigenvalues for 
several plane and three-dimensional self-conjugate equations — 
the linear system is changed to be plane or three-dimensional 
(for example, finding the eigen functions of the equation 
etc. 

* 

Note . Let us assume we are solving an experimental problem 
of finding the eigenvalues and the eigenvectors for the equation 

Cy is the n-dimensional vector; A - square matrix). The solution 
of this problem may be modeled very simply if the diagonal terms 
have one sign in the matrix A. This can be achieved by adding the 
corresntriding positive constant C to the diagonal terms. In addi- 
tion, the eigenvectors are not changed, and the eigenvalues in- 
crease by C [our equation is replaced by the equation (A + C)y * 

(), +- C) V J. 

We give this in the problem examined. 

The answers obtained wil3 contain inaccuracies following both 
from the unreliable estimate of the measurements, and from the dis- 
tortion which we assume, as well as by the replacing of the differ- 
ential operator by the finite-difference operator. For a plane 
problem to be solved by the grid method, we may replace an arbi- 
trary contour by a polygon contour. However, we may greatly improve 
the results by regarding the distortions we introduced as perturba- 
tions of the operator and by determining the influence of these 
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perturbations on the result obtained. This pertains not only to 
the method given, but to any other experimental or numerical method 
of solving these problems. 


Section 6 


Let L be a self-conjugate operator and we shall approximately 
solve the operator equation 


Lu — >u — 0 

under the normalization condition 


(u,u)~l. 


(23) 


(24) 


Actually, we replace the operator L by the approximating opera- 
tor (we replace the differences by the derivative in the differ- 
ential operator, the integrals by the sums in the integral operator, 
etc.), where (in the previous example 5 L is the difference 

between the differential and the finite difference operators, which 
is expressed by means of the difference of higher orders). 


We have found the approximate eigenvalue ^ and the approxi- / 6 2 

mate normalized eigen element u 1 

L i —L 4- IL, Uj ~ u -+* &u, 
which satisfy the equation 

L, u, — - u -= e (25) 

(the nonzero right side is determined by calculation errors, mea- 
surement errors, etc.). 

Substituting (23) from (25) and disregarding terms of the 
second order of magnitude with respect to 5 A and 5 u, i.e., replac- 
ing the difference of the first variation, we obtain 
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Ukt iLu — u — hu 


i. 


(26) 


— ' — i 


Multiplying both sides of C 26 ) by u and using the equations 
(23) and (24), and in addition, due to the self-conjugate nature ! 

of L, we obtain 1 

1 

{Lu — Au) + (iLu, u) — h.= («, u), 1 

" ii 

' i 

A 

from which we have I 


* = (''Lu, u) — (i, «). 


(27) 


The first term on the right side is the regular expression 
for the perturbation of the eigenvalue. The second is determined 
by the inaccuracies of solving the "perturbed” equation. Finding 
SA and designating by yj — dA, we obtain a more precise value 
of A. Equation (.26), which is solved with a relatively small 
accuracy, determines the correction d . 

In problems for finding the ei-envalues of the differential 
operators, the boundary values are included ~ the zero values 
at the boundary. Let us use the symbol & 2 t0 designate the varia- 
tion caused by varying the boundary. For an approximate solution, 
if we assume the distortion of the boundary (for example, in the 
plane problem we replace the curvilinear plane contour by a con- 
tour with several corners), then the termdgA appears in the right 
side of C27) — the variation of the eigenvalue caused by deforma- 
tion of the boundary. 

R. Courant developed (see, for example, Part 1 of "Methods 
of Mathematical Physics") a variational theory of eigenvalues 
based on the method of the minimum-maximum. However, this theory 
gives only a qualitative estimate — the increase or decrease of 
the eigenvalues, i.e., only the sign of wh *l e the regular 
elementary methods of the variational value give the values 5 
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i.e., they determine (within an accuracy of the higher variations) 
the numerical value of the distortion Introduced by the small 
deformation of the boundary, and make it possible to refine the 
result. 

Let us assume L is the self-conjugate differential operator 
with certain zero boundary data. 

Let us use /(«)=(£«, u), to designate 4(u)=S(u)~X(u,u) . Equa- 
tion (23) is the Euler equation 3/(«)s=0s Under the isoperi- 
metric condition (u,u)» 1 or <V, i«) 0 . The normalized eigen /6 3 

function is the extremum of the equation and the eigen- 

value J («*.“)• 

We shall vary the boundary. Then V'— V- . Since the normal- 

ization condition retains its form, 

% t (u,n)~0, t . e. 

The last side of this equation is calculated on the basis of 
the regular rules for the variation of the functional from the ex- 
tremum with a variation of the boundary. 

Example 1 . L is the Sturm-Liouville operator. We have the 
equation (17) under the condition u(a) ■ u(b) ■ 0 and the normali- 
zation 

ft 

( u,u t = ) tr dx — 1. 

Here 


4 

J (u) — (Lu,u) - I (Ru * Pll)dx, 

t 

y,(u)= f (Ru * 

Setting F Ru* •• Pu : — ‘>u , we have 
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F.> = 2Ru, F— u’F". — — Ru* i- /»«--> u-. 

t 

If the finite abscissas a and b obtain the variation *a and Ab, 
*-(/-**' '*> > ■-(F-U'F,,;**. 

At the point square x ■ a, x ■ b, we have 

F — u F„> — — Ru' 1 , 

Therefore 

*t >• « *<«) la’ (a)f — R (b) [«’ (6)}- *b. (28) 

Since 


R (a) > 0, R b)> 0, 

then at Jt/?>Q,&6<0 [i.e., with compression of the Integral 

(a,b)] 

^>C, 

i.e., the eigenvalues Increase. Thus, the qualitative side of 
the variational series of eigenvalues is completely obtained. 

Formula (28) gives the quantitative estimate. 

Example 2. Let us assume L ■ V is the plane Laplace opera- 
tor applied to the functions u(x,y) given in the plane domain Q /6fr 

with zero boundary values at the boundary q of the domain Q. Then 

•'(«)= I j (v [({*)]■ 

y > w - J J L( S )’ ■*• (t)’ “ J * 

Q 

Let us vary the boundary q In the vicinity of its point A, and we 
assume ^ is the variation of the area Q. Then 

i,/,(u)=- [grad u (4)]* 
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and, consequently, 

a,X = — (graduM)! 3 *^. 

Thus, the qualitative estimate of fi 2 * 18 obtained: at ^i<°. *i*>0. 

i.e., the eigenvalues of the operator V increase, if the domain Q 
is reduced, changing into the eigen part. However, formula (28) 
gives the quantitative estimate. 

Let us assume we have found the eigenvalue for A and the 
eigen functions of u for the domain which approximates the 
domain Q (for example, for the domain consisting of the quadrile 
squares, if we solve experimentally the problem for the plane 
square grid). Let us divide the band between Q and into small 

parts of the area * = L 2 » which are adjacent to the 

point Aj, i ■ 1, 2, . . . , n of the boundary q. With the transition 
from to Q, the eigenvalue \ receives the increment 5 2 A which ap- 
proximately equals 


■' ^ [irad (!(.■? )] ^ , 
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